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We report on a search for Standard Model f-channel and s-channel single top quark production 
in pp collisions at a center of mass energy of 1.96 TeV. We use a data sample corresponding to 
162 pb -1 recorded by the upgraded Collider Detector at Fermilab. We find no significant evidence 
for electroweak top quark production and set upper limits at the 95% confidence level on the 
production cross section, consistent with the Standard Model: 10.1 pb for the t-channel, 13.6 pb for 
the s-channel and 17.8 pb for the combined cross section of t- and s-channel. 

PACS numbers: 14.65.Ha, 12.15.Ji, 13.85.Rm 



In pp collisions at 1.96 TeV, top quarks are predomi- 
nantly produced in pairs via strong interaction processes. 
Within the Standard Model (SM), top quarks are also 
expected to be produced singly by the electroweak in- 
teraction involving a Wtb vertex At the Tevatron, 
the two relevant production modes are the t- and the 
s-channel exchange of a virtual W boson. The mea- 



surement of the single top cross section is particularly 
interesting because the production cross section is pro- 
portional to |Vtb| 2 , where Vtb is the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element which relates top and 
bottom quarks. Assuming three quark generations, the 
unitarrty of the CKM matrix implies that Vtb is close to 
unity Q|. The most recent next-to-leading order (NLO) 
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calculations, assuming |Vt&| = 1, predict cross sections of 
(1.98 ±0.25) pb for the f-channel and (0.88 ±0.11) pb for 
the s-channel mode at ^fs = 1.96 TeV Q. Using these 
predictions, a measurement of the single top cross section 
will allow for a direct determination of \Vtt\- Single top 
searches test also exotic models which predict anoma- 
lously altered single top production rates [4|. Moreover, 
single top quark processes result in the same final state 
as the Standard Model Higgs boson process WH — > Ivbb 
and therefore impact future searches for the Higgs boson 
at the Tevatron 5]. In this article we report results of 
the first search for single top production in Run 2 at the 
Tevatron. Results of searches for single top production 
at = 1-8 TeV (Run 1) can be found in Refs. 0,0. 

The experimental signature of single top events con- 
sists of the W decay products plus two or three jets, in- 
cluding one b quark jet from the decay of the top quark. 
To suppress QCD multijet background we select only 
W — ► fw^ and W — > ev e candidates. In s-channel events 
we expect a second b quark jet from the Wtb vertex. In 
t-channel events a second jet originates from the recoil- 
ing light quark and a third jet is produced through the 
splitting of the initial-state gluon into a bb pair. Mostly, 
this third jet escapes detection, since it is produced in the 
high pseudorapidity (77) regions and at low transverse en- 
ergy (E T ) 0. 

This article describes two analyses: (1) a combined 
search for the t- plus s-channel single top signal, (2) a 
separate search, where we measure the rates for the two 
single top processes individually. The data sample corre- 
sponds to an integrated luminosity of (162 ± 10) pb -1 
collected with the upgraded Collider Detector at Fer- 
milab (CDF II), which is described elsewhere 0- The 
common event preselection of our two analyses resem- 
bles closely the one used in the CDF measurement of the 
it cross section reported in Ref . 0| ■ We accept events 
with evidence for a leptonic W decay: (a) missing trans- 
verse energy ]£t > 20 GeV from the neutrino and (b) 
an isolated central electron with Et > 20 GeV or an 
isolated central muon candidate with pt > 20 GeV/c. 
An electron or muon candidate is considered isolated if 
the non-lepton Et in an rj-cj> cone of radius 0.4 centered 
around the lepton is less than 10% of the lepton Et or 
Pt- To remove dilepton events from tt-production and 
leptonic Z boson decays, we accept events with only one 
well identified lepton. In addition, we veto events if we 
find a second, loosely identified lepton candidate that 
forms an invariant mass with the primary lepton between 
76 GeV/c 2 < M u < 106 GeV/c 2 . The jet reconstruction 
uses a fixed cone of radius AR = \J A(j) 2 + Arf = 0.4. 
We count jets with transverse energy Et > 15 GeV and 
\rj\ < 2.8. We only accept W + 2 jets events. At least 
one of these jets must be identified as likely to originate 
from a b quark (6-tag) by requiring a displaced secondary 
vertex within the jet as measured using silicon tracker in- 
formation. The effective coverage of the 6-tagging ranges 



TABLE I: Event detection efficiencies in %. 



Process 


Combined 


1-tag 


2-tag 


t-channel 


0.89±0.07 


0.86±0.07 


0.007±0.002 


s-channel 


1.06±0.08 


0.78±0.06 


0.23±0.02 



up to |?7| < 1.4. 

To optimize our sensitivity, we apply a cut on the in- 
variant mass Mi u b of the charged lepton, the neutrino 
and the 6-tagged jet: 140 GeV/c 2 < M lvb < 210 GeV/c 2 . 
The transverse momentum of the neutrino is set equal to 
the missing transverse energy vector i? T ; p z {v) is ob- 
tained up to a two-fold ambiguity from the constraint 
Mg v = Mw- From the two solutions we pick the one with 
lower |p z (z/)|. If the p z {v) solution has non-zero imagi- 
nary part as a consequence of resolution effects in mea- 
suring jet energies, we use only the real part of p z [y). For 
the separate search, we subdivide the sample into events 
with exactly one &-tagged jet or exactly two 6-tagged jets. 
For the 1-tag sample, we require at least one jet to have 
Et > 30 GeV. We determine the total event detection 
efficiency e cv t for the signal from events generated by the 
matrix element event generator MadEvent followed 
by parton showering with PYTHIA [T3| and a full CDF 
II detector simulation. MadEvent features the correct 
spin polarization of the top quark and its decay prod- 
ucts. For i-channel single top production we generated 
two samples, one b + q — ► t + q' and one g + q — > t + b + q' 
which we merged together to reproduce the pt spectrum 
of the b as expected from NLO differential cross section 
calculations. This is an improved model compared to the 
Pythia modelling used in the Run I analyses. The event 
detection efficiency e ov t includes the kinematic and fidu- 
cial acceptance, branching ratios, lepton and 6-jet identi- 
fication as well as trigger efficiencies. We combine e cv t as 
given in Table U with the cross sections predicted by the- 
ory [3j and thereby obtain the number of expected single 
top events listed in Table ITT1 

We distinguish between two background components: 
it and non-top background. We estimate the it back- 
ground based on events generated with PYTHIA, nor- 
malized to the theoretically predicted cross section of 
a{ti) = 6.7I.o;9 pb 13 . The primary source (62%) of 
the non-top background is the W+heavy flavor processes 
qq' — > Wg with g — > bb or g — > cc, and gq — > Wc. Ad- 
ditional sources are "mistags" (25%), in which a light- 
quark jet is erroneously identified as heavy flavor, "non- 
W"(10%), e.g. direct bb production, and diboson (WW, 
WZ, ZZ) production (3%). The non-W and mistag frac- 
tions are estimated using CDF II data. The Fy±heavy 
flavor rates are extracted from ALPGEN 0] Monte 
Carlo (MC) events normalized to data [llij . The diboson 
rates are estimated from PYTHIA events normalized to 
theory predictions |l5j . The numbers of expected sig- 
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TABLE II: Expected number of signal and background events 
compared with observations. 





(~!mnbinpd 


l Lag 


2-tag 


t-channel 
s-channel 


2.8 ±0.5 
1.5 ±0.2 


2.7 ±0.4 
1.1 ±0.2 


0.02 ±0.01 
0.32 ±0.05 


ti 

non-top 


3.8 ±0.9 
30.0 ±5.8 


3.2 ±0.7 
23.3 ±4.6 


0.60 ±0.14 
2.59 ±0.71 


Total Background 


33.8 ±5.9 


26.5 ±4.7 


3.19 ±0.72 


Total Expected 


38.1 ±5.9 


30.3 ±4.7 


3.53 ±0.72 


Observed 


42 


33 


6 
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FIG. 1: Q -r] distributions for a) MC templates normalized to 
unit area, b) Data in the 1-tag sample (33 events) vs. stacked 
MC templates normalized to the SM prediction. 



nal and background events are summarized in Tabic ITT1 
Having applied all selection cuts we observe 42 events for 
the combined search, 33 events in the 1-tag sample and 
6 events in the 2-tag sample. Within the uncertainties, 
the observations are in good agreement with predictions. 

To extract the signal content in data, we use a maxi- 
mum likelihood technique. We separate t- and s-channel 
events by using the Q ■ r\ distribution which exhibits a 
distinct asymmetry for t-channel events, see Fig. QJi. Q 
is the charge of the lepton and r\ is the pseudorapidity 
of the untagged jet. In Fig. we show the data versus 
stacked MC templates weighted by the expected number 
of events in the 1-tag sample. The separate search de- 
fines a joint likelihood function for the Q ■ i] distribution 



in the 1-tag sample and for the number of events in the 
2-tag sample. 



£ s i g ((7i, ... , 0-4; Si, ■■■ ,87) 



-Hi . ..n d ^bin Mfc n k 



rid 



fc=l 



i>k 



3 = 1 i=l 

Four processes are considered and labeled by the index 
j: t-channel (j = 1), s-channel (j = 2), ti (j — 3), and 
non-top ( j = 4) . The corresponding cross sections are de- 
noted Oj . The background cross sections are constrained 
to their SM prediction ctsmj by Gaussian priors of width 
A 3 = 23%ctsm,3 for ti and A4 = 20%ctsm,4 for non-top. 
The index "sig" denotes the signal process, which is s- or 
t-channel, respectively. The fik are the mean number of 
events in bin k of the Q ■ 77 histogram (Nbi n = number 
of bins), while fid is the mean number of events in the 
2-tag sample, and nd are the event numbers observed 
in data, respectively. Seven sources of systematic un- 
certainties are considered in the likelihood function: (1) 
jet energy scale (JES), (2) initial state radiation (ISR), 
(3) final state radiation (FSR), (4) parton distribution 
functions (PDF), (5) the choice of signal MC generator, 
(6) the top quark mass, (7) trigger, identification and 
^-tagging efficiencies and the luminosity. The relative 
strength of a systematic effect due to source i is parame- 
terized by the variable Si. Systematic effects change the 
acceptance and influence the shape of the Q ■ r\ distri- 
bution. When calculating f^k/d we take the systematic 
shifts in the acceptance and in the shape of the tem- 
plate histograms, and their full correlation into account. 
All variables except the signal cross section a S i g are con- 
strained to their expected values by Gaussian functions 
G(x;xq, A x ) of mean Xq and width A x . The largest un- 
certainties are on the 6-tagging efficiency (7%), luminos- 
ity (6%), top quark mass (4%) and JES (4%). The effect 
of uncertainty in the JES is evaluated by applying energy 
corrections that describe ±1 a variations. Systematic un- 
certainties due to the modeling of ISR and FSR are ob- 
tained from MC samples that describe variations in these 
effects. To evaluate the uncertainty associated with the 
choice of a specific parametrization of PDF we investi- 
gated several PDF sets and took the maximum deviation 
(MRST72) from our standard PDF set (CTEQ5L). We 
estimate the uncertainty associated with the choice of 
single top MC generator using samples generated with 
TopReX 0| . The values of acceptance uncertainties for 
the single top processes are summarized in Table Hill 

To measure the combined t-channel plus s-channel sig- 
nal in data, we use a kinematic variable whose distri- 
bution is very similar for the two single top processes, 
but is different for background processes: Ht, which is 
the scalar sum of Ifix and the transverse energies of the 
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TABLE III: Fractional changes in e cv t of single top processes 
in %. e t rig is the trigger efficiency, eiD the lepton identification 
efficiency. 



i 


Source 


f-channel 


s-channel 


Combined 


1 


JES 


+2.4 

-6.7 


+0.4 
-3.1 


+ 0.1 

-4.3 


2 


ISR 


±1.0 


±0.6 


±1.0 


3 


FSR 


±2.2 


±5.3 


±2.6 


4 


PDF 


±4.4 


±2.5 


±3.8 


5 


Generator 


±5 


±2 


±3 


6 


Top quark mass 


+0.7 
-6.9 


-2.3 


-4.4 


7 


Etrigj £ id, luminosity 


±9.8 


±9.8 


±9.8 




100 150 200 250 300 350 400 450 
H T [GeV] 

FIG. 2: Ht distribution for data (42 events) in the combined 
search compared with smoothed MC predictions for signal and 
background. 



lepton and all jets in the event. We use a likelihood func- 
tion similar to that in the separate search. One difference 
is that in the combined search the Ht distributions are 
smoothed. In Fig. [3] we show the Ht distribution ob- 
served in data compared to the SM prediction. 

We perform Monte Carlo experiments to estimate our 
a priori sensitivity assuming the SM signal cross sections. 
For each experiment we integrate out all nuisance param- 
eters (i.e., all variables except o~ S ig) from the full likeli- 
hood function and thereby construct the marginalized 
likelihood C* ig {a S i g ). C* ig is normalized and interpreted 
as posterior probability density function p{a S i g ). We cal- 
culate the upper limit at the 95% C.L. using a Bayesian 
method assuming a prior probability density, which is 
if a S i g < and 1 if o~ S i g > 0. The median of the expected 
upper limits defines our sensitivity and is stated in Ta- 
ble El We calculate the posterior probability densities 
for CDF II data and obtain the most probable values 
(MPV) and highest posterior density (HPD) intervals Q 
as given in Table llVl Within the statistical uncertainty 
these results are compatible with SM predictions. We 



TABLE IV: Upper limits at the 95% C.L. and most probable 
values (MPV) of single top cross sections in pb. 





t-channel 


s-channel 


Combined 


expected limit 
observed limit 
MPV ± HPD 


11.2 
10.1 

0+ 47 


12.1 
13.6 
4 fi+ 38 

4 - D -3.8 


13.6 
17.8 

7 7+5.1 
'•'-4.9 



find upper limits of 10.1 pb at the 95% C.L. for the t- 
channel cross section and 13.6 pb for the s-channel. For 
the combined search we find an upper limit of 17.8 pb at 
the 95% C.L.. 

In summary, we find no significant evidence for elec- 
troweak single top quark production in (162 ± 10) pb -1 
of integrated luminosity recorded with CDF II. We set 
the first limits on single top cross sections in pp collisions 
at ^/s — 1.96 TeV in Run 2 at the Tevatron. If com- 
pared with Run 1 results the upper limits on ^-channel 
and s-channel single top quark production are consider- 
ably improved by 28% (i-channel) or 20% (s-channel), 
respectively. We have introduced improved Monte Carlo 
modeling for single top and a fully Bayesian treatment of 
systematic uncertainties in the likelihood function which 
are important steps for future analyses aiming for the 
observation of single top quark production. 
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